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Relations between photon scattering, entanglement and multi-mode detection are investigated. 
We first establish a general framework in which one- and two-photon elastic scattering processes 
can be discussed, then we focus on the study of the intrinsic entanglement degradation caused by a 
multi-mode detection. We show that any multi-mode scattered state cannot maximally violate the 
Bell-CHSH inequality because of the momentum spread. The results presented here have general 
validity and can be applied to both deterministic and random scattering processes. 
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I. INTRODUCTION 

During the last fifteen years, the concept of entan- 
glement has passed from the inhospitable realm of the 
"problems in the foundation of quantum mechanics" [J, 
to the fashionable fields of quantum information and 
computation 2\ . Experimental availability 3] of reliable 
sources of entangled quantum bits (qubits) in the optical 
domain has been one of the major boosts for this rapid 
transition. Polarization-entangled photon pairs produced 
in a spontaneous parametric down conversion (SPDC) 
process Q have thus become an essential tool in experi- 
mental quantum information science . 

Robustness of photon polarization-entanglement is im- 
portant when propagation over long distances of entan- 
gled states is required. However, during propagation in 
a medium photons may be depolarized by scattering pro- 
cesses due to imperfections. Further, the recently ob- 
served transfer of quantum entanglement from photons 
to a linear medium and back to photons 0, has also been 
described as an elastic scattering process |3- In general, 
polarization-dependent scattering may affect more or less 
the entanglement of a photon pair depending on the prop- 
erties of the scattering medium |^ ■ 

In this paper we show that independently from the 
details of the scattering process, the polarization entan- 
glement of a photon pair is unavoidably degraded by a 
multi-mode measurement |9|. Moreover, we show that, 
contrary to common belief, it is not always possible to 
build a 2 X 2 reduced density matrix (or 4 x 4 in the case 
of entangled pairs) to describe the scattered state. This 
fact was recently recognized by Peres and coworkers [Tol | 
for the case of a single polarization qubit, as opposed to 
the case of two-photon entangled qubits considered here. 
They circumvented this difficulty by introducing, at a 
certain stage of their calculations, the unphysical con- 
cept of longitudinal photons. In this way they were able 
to build a 3 X 3 effective reduced density matrix whose 
physical content can be obtained from a family of positive 
operator- valued measures (POVMs [Tl|). 

In this paper we follow a somewhat different approach. 
We show that the difficulties one encounters in trying 
to define a reduced density matrix, are also due to the 



troubling determination of the effective dimension of the 
Hilbert space in which the two-photon entangled state 
can be represented. In fact, due to the transverse nature 
of the free electromagnetic field, in a multi-mode scat- 
tering process the momentum and the polarization de- 
grees of freedom become entangled. Therefore, after the 
scattering, the two-photon polarization-entangled state 
is no longer confined to a 4-dimensional Hilbert space, 
but spans an higher-dimensionality space. We show that 
under certain circumstances this dimension may remain 
bigger than 4 even after tracing out the momentum de- 
grees of freedom. In this case a traditional Bell-violation 
measurement setup may fail to reveal a maximal viola- 
tion of the Bell-CHSH inequality. However, the detailed 
investigation of the entanglement properties of this kind 
of high-dimensional scatteredphoton-pair 0, ITsL Q| , is 
left to a forthcoming paper 
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This paper is structured as follows. In Sec. II we de- 
fine the general framework for the discussion of one- and 
two-photon scattering processes. In Sec. HI, starting 
from an experimentalist point of view, we adopt an op- 
erational way to define the photon polarization states in 
terms of the physically detected states. Furthermore, we 
calculate the classical Jones matrix corresponding to 
an arbitrarily oriented polarizer and use it to introduce 
a suitable field creation operator to build explicitly the 
physical photon states. In this way we can introduce a 
2x2 effective density matrix by a projection operation 
on the one-photon physical states. The concept of effec- 
tive density matrix is extended to polarization-entangled 
photon-pair states, in the first part of Sec. IV. Then, 
in the second part of Sec. IV we quantify the entangle- 
ment degradation by calculating the effect of a polariza- 
tion detection mismatch on the mean value (.Bchsh) of 
the Bell-CHSH operator 17] for the set of the Schmidt 
states |l8j . Finally, we summarize our main conclusions 
in Sec. V. 
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II. ONE-PHOTON SCATTERING 
A. One-photon states and their representation 

To understand the effects that a scattering process 
may have on a single photon wave packet, we start by 
introducing a proper notation. Let A{r,t) denote the 
transverse part of the electromagnetic potential. In the 
Coulomb gauge the quantization procedure can be car- 
ried out straightforwardly to obtain the correspond- 
ing field operators A{r,t). At any given time t, A{r,t) 
may be expanded in terms of the Fourier series 
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(1) 



where to = |k|, ak(i) ■ k = 0, and natural units (c ~ h ~ 
1) have been used. With H. we denote the quantization 
volume. For any given k we may define a set of three real 
orthogonal unit vectors k = k/|k|, e'^\ e^^^ such that 



.(1) a(2) 



C^=0, eWx4^)=k 



(2) 



It is convenient to introduce the creation operators 
of a photon with momentum k and helicity s = ±1 as 



where the complex vectors f]^"' are defined as 



(3) 



f(«) 



v/2 



,(2: 



" V2 

They satisfy the canonical commutation rules 



(. = ±1). (4) 



aks{t),a'l,^,{t) = (5kk"5, 



(5) 



Then, an one-photon helicity basis state can be written 
as 

|k,s(k)) ^aL|0), (6) 

where |0) denotes the vacuum state of the free elec- 
tromagnetic field. This state represents a photon with 
momentum k and helicity s(k) = ±1, that is |k, s(k)) 
is an eigenstate of the momentum operator K = 
Sk' s' k'aj^,^,ak's' with eigenvalue k: 



K|k,s(k)> =k|k,s(k)). 



(7) 



Because of Eq. O, the onc-photon-space basis states 
|k, s(k)) are often said to be the direct products of mo- 
mentum and polarization states and this is usually stated 
by writing 



|k,s(k))^|k)®|s(k)). 



(8) 



However, one should realize that the left side of this equa- 
tion is a genuine QED state, while the right side is just 



one of its possible representations. To be more specifics 
in QED there is no such a thing like a momentum cre- 
ation operator that creates from the vacuum a momen- 
tum state |k), nor is there a polarization creation opera- 
tor which create a polarization state |s(k)). Nevertheless 
the "momentum (8) polarization" representation for the 
one-photon states has several advantages. For this rea- 
son we now present formally a few elementary facts about 
this representation. 



1. Some elementary facts 

Let us consider two one-photon states and |0q) 
which are eigenstates of the momentum operator with 
eigenvalues k and q respectively: 



|7/^k> - V^+(k)|k,l)+7^^(k)|k,-l), 

\4>^) = <^+(q)|q,l) +0_(q)|q,-l). 
They can be represented as 



l^k 



(9) 



(10) 



where the symbol stands for "is represented by" , and 
we have used the subscripts k and q for the representa- 
tions of iT/^k) and |0q) respectively, to stress the fact that 
they belong to different vector spaces. Since 



we can write 

IV^k) =^+|+)k + V'-|-)k, 

where we have defined 



^k = 



jk = 



(11) 



(12) 



(13) 



In the equations above we have used the parenthesis sym- 
bols |±)k for the basis vectors instead of the usual bracket 
|-) ones to underline the fact that they are not genuine 
QED states, that is they are not created from the vac- 
uum by any bosonic creation operator. Then, from now 
on, we shall reserve the bracket symbols |-) only for the 
truly QED states. 

The two vectors displayed in Eq. (|13() form a basis 
for a two-dimensional "polarization" space Tik- Ex- 
actly in the same way we can introduce a basis Sq = 
{|+)q, | — )q} for the space Tiq relative to the (arbitrary) 
momentum q. More generally, let if be a iV-dimensional 
set of momenta: K — {ki, . . . , kAr}. Then we can 
build a 2iV-dimensional space TL spanned by the basis 
S = Ski U • • • U 6k„ as 



'H = Hki 



>Wk. 



(14) 
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It is clearly possible to introduce other bases in H.. For 
example let {|ki)}, (i = 1,...,A^) and {|+),|— )} be a 
iV-dimensional and a 2-dimensional sets of basis vectors 
respectively, represented by 



|k.) = 



(15) 



where all elements of |ki) are zero with the exception of 
element i which is equal to 1. Then the set given by the 
vectors 



= l,...,iV, 

s = ±1, 



(16) 



is an orthonormal and complete basis of 7i. Now, the 
factorization "momentum ® polarization" in Eq. H16() is 
absolutely legitimate without introducing any "ad hoc" 
creation operator. We would like to stress once again 
the fact that the genuine QED states |k, s(k)) are not 
factored, just their representation |k, s). Moreover, the 
helicity label s in Eq. (|16l) does not depend on the mo- 
mentum k, so it is not longer a secondary variable |lClj| . 
Note that we have chosen to work with the helicity eigen- 
states |k, s(k)) but we could have worked equally well 
with the linear polarization basis vectors |k, Q;(k)), where 
a = x,y. 



2. Tracing over the momenta 

The reason for which we have introduced the factoriz- 
able basis vectors |k, s) in the previous subsection, is that 
we want to be able to separate the momentum degrees 
of freedom from the polarization ones. This separation 
becomes important when the photons are subjected to 
polarization tests only, and the momenta become irrele- 
vant degrees of freedom. Then, it is a standard practice 
[ll| to introduce a reduced density matrix obtained from 
the complete one, by tracing over some subset of the 
whole set of the field momenta. 

Let p be the density matrix describing the state of a 
single photon. We can represent /5 in 7i by defining the 
iV^ matrices R{n,m) (2 x 2) whose elements are calcu- 
lated with the usual QED rules: 

Rss'{n,m) = Zi(k„,s(k„)|p|k™,s'(k„)), (17) 

where Zi is an arbitrary normalization constant and 
(s,s' = ±1). These matrices are the building blocks of 
the total representation of p: 



( 



P = 



R++{1,1) 
R-+{1,1) 



i?+_(l,l) 
R—{IA) 



R{1,N) 



R{N, 1) 



(18) 



R{N,N) J 



From now on, we choose Zi such that Trp = 1. Now, if we 
want to calculate the redwced density matrix pf^, obtained 
by tracing out the momentum degrees of freedom, what 
we have to do is simply to take all the diagonal blocks 
R{n, n) and sum them together: 



Pf 



N 



R{n, n) 



pf- 



(19) 



The matrix pf is a well defined 2x2 matrix and there 
is no ambiguity in its determination. We have used the 
subscript 1 in writing pf to distinguish it from the gen- 
eralized reduced density matrix p^ we shall introduce 
in the next subsection. However, it was recently shown 
that a reduced density matrix introduced as above, has 
no meaning because it does not have definite transfor- 
mation properties Besides this fact, there are more 
practical reasons for which pf (at least a single one) can- 
not be introduced, as we shall see in the next subsection. 
For the moment let us see very shortly how the tradi- 
tional approach works and why it may fail. If one deals 
with operators that are factorizable in the basis {|ki, s)}, 
that is if they can be represented as 



N 



1- |X)(X|®^|k„)(k„|, 



(20) 



where the right side of this equation is, in fact, the direct 
product of a 2 X 2 matrix (polarization part, momen- 
tum independent), times a, N x N matrix (momentum 
part, polarization independent). Then the average value 
{X) = Ti{pX) can be certainly written as 



{X)^Tr{pf\X)iX\) 



(21) 



and, apparently, there are no problems in using pf. The 
key point here, is that the operator representation with 
respect to the factorizable basis {|k„,s)} may have no 
meaning at all because the polarization state |X) in Eq. 
(|20|l docs not depend on k while the polarization state 
of a photon is always referred to a given momentum, as 
shown in Eq. (j3Jl. Then, in a multi-mode process where 
many values of the photon momentum are involved, it be- 
comes impossible to define a unique 2x2 reduced density 
matrix (|X)(X| in Eq. if^ ). instead one must deal with 
a different matrix for each different momentum. This will 
be explicitly shown in the next subsection. 



B. One-Photon Scattering 

Let us consider now a one-photon state approxima- 
tively represented by the monochromatic plane wave 
|ko, so(ko))in, and suppose that it is elastically scattered 
by a certain medium. The final state {ipf) of the photon 
after the scattering process can be written in terms of 
the output states |k, s(k)) as 



i^/) = E E ^^(k)|k,.(k)), 



(22) 



keK s=±l 
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where V's(k) = (k, s(k)|ko, so(ko))in represent the prob- 
abihty ampHtude that the photon is scattered in a state 
with momentum k and hehcity s(k), where |k| = ujq. 
With K we have denoted the set of ah scattered modes. 
By inspecting Eq. (|22l) one immediately reahzes that 
in the scattered state {ipf) the polarization and momen- 
tum degrees of freedom are entangled because the helicity 
states |s(k)) depend explicitly from k and the two sums 
are not independent. In other words, in general, {ipf) is 
not an eigenstate of the linear momentum operator K. 

Now suppose that we measure some polarization prop- 
erty of the scattered photon, regardless of its momen- 
tum. To be more specific let us assume that a polar- 
ization analyzer (a dichroic sheet polarizer or a crystal 
prism polarizer) is present in a plane perpendicular to 
the wave-vector kg of the impinging photon and that we 
collect, with a photo-detector, all the light coming from 
the scattering target within a certain angular aperture 
0£). With this experimental configuration we test only 
the photon polarization, irrespective of the plane-wave 
mode k where the photon is to be found within Qo- We 
call this kind of experimental arrangement, a multi-mode 
detection scheme. 

The act of measurement can be quantified by calcu- 
lating, for example, the mean value of the operator P 
defined as: 

P= ^ ^P,,,(k)|k,.(k))(k,.'(k)|, (23) 



where 



p(k) 



p++(k) p+_(k) 

F_+(k) P__(k) 



(24) 



is a 2 X 2 Hermitian matrix such that — P and Ko is 
the set of the detected modes: 

{\^eKD\ Ik-kol/wg < IcosGdI}. 

Moreover we assume, for definiteness, dim K£i = N and 
Kd C K. It is easy to see that P is a projector (P^ — P), 
and that it commutes with the momentum operator K: 



[P,K] =0. 



(25) 



It is also easy to understand that the equation above is 
the equivalent, in the QED context, of the factorizability 
condition expressed in Eq. 120() . The matrix elements of 
P are 



Pss'{n,m) = (k„,s(k„)|P|k„,s'(k™)) 



(26) 



From Eq. H2()|l is clear that P has a block-diagonal shape: 



P = 



P++(l) P+_(l) 
P-+(l) P-{1) 



(27) 



PiN)J 



so that only the corresponding diagonal blocks of p will 
enter in the calculation of (P) . Explicitly 



N 



(P) =^Tr(P(n)P(n)), 



(28) 



where Eq. H18|l has been used and R{n) = R{n, n) for 
short. Then, the form of the operator P in Eq. (|27|l nat- 
urally leads to the introduction of a generalized reduced 
density matrix {2N x 2A^) defined as: 



P(l) 



p« =diag{P(l),...,P(7V)} = 



(29) 

Thus, differently from Eq. H19|l were we have found a 
unique 2x2 matrix, we have now to deal with a set of 
N different 2x2 matrices {P(l), . . . , R{N)}: one matrix 
R{n) per each mode k„ of the field. This is due to the 
form of the operator P. Only in the hypothetical case in 
which the matrix elements Pss' (k) of P were independent 
from k, we would obtain again a unique 2x2 density 
matrix defined as pf = J2n=i ^("-)- 

What does all this mean? We can get some insights 
into the above equations, by rewriting them in the basis 
B — ;Bki U • • • U ;Bk„ introduced previously. Let us start 
by defining the normalized state |(/))k in the basis as 



VC(k) ,=±1 



(30) 



where the complex amplitudes V's(k) are defined in Eq. 
(ESl and 



C(k)^|V+(k)p + |^_(k)p. 



(31) 



Then we can represent the one-photon state \i})f) in Eq. 
()22|l . restricted to Kd, as 



(32) 



ke-fsTr 



Then, by combining Eq. (flTjl with Eq. |23) we find 

Rss'{n,m) = ZiV's(k„)V'*,(k„). (33) 

From the equation above and the normalization condition 
TrP = 1, it follows that 



(34) 



keA'i 



It should be clear now that if we define the pure-state 
density matrix p^Ck) as 



/50(k) ^ 



k kl' 



(35) 
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then we can write the generalized reduced density matrix 
as the direct sum, as opposed to the ordinary sum in 
Eq. l(TT?]l . of all the sub-matrices ^^(k): 



= |,/.)k«;(k)k(</'|, 

were the statistical weight function w(k) is equal to 

C(k) 



(36) 



wik) = ZiC(k) 



(37) 



and 



E 



w(k) = 1. So we find that the price to pay 



for introducing a polarization basis B, is that we have to 
deal with a generalized reduced density matrix which is 
expressed as a direct sum instead of an ordinary sum. Of 
course, this fact is entirely due to the incomplete nature 
of the polarization representation; there are no problems 
in writing as an ordinary sum in the QED formalism. 
In such a context we have to use the complete momentum 
eigenstate basis {|k, s(k))} in order to write 



P 



«;(k)|k,0(k))(k,(/.(k)|, 



(38) 



keKr 



where we have defined the momentum eigenstates 
|k,0(k)) as 

|k, 0(k)) s -1= J2 V'^WIk, s(k)). (39) 

Not surprisingly, has the same "momentum-diagonal" 
structure as P in Eq. H23|l . It clearly represents a mixed 
state because (p^)^ 7^ . Eq. (jSHl can be interpreted by 
saying that from the observer point of view (who suppos- 
edly cannot measure the photon momentum), things go 
as if the scatterer were a thermal source emitting photons 
in the states |k, s(k)) with probabilities ti;(k). 

Thus, we have shown that it is impossible to extract 
from p^ in a straightforward way a unique 2x2 reduced 
density matrix; rather we have obtained a iV-dimensional 
set {R{1), . . . , R{N)} of 2 X 2 matrices, each of them 
well defined in its proper Hilbert space specified by the 
momentum k. This is a direct consequence of the fact 
that helicity of a photon can be defined only with respect 
to a given momentum It is worth to note that this 
result is not in contrast with existence of a 4 x 4 density 
matrix for two polarization-entangled photons emitted 
in a SPDC process. In fact, in that case, each photon is 
emitted in one mode of the electromagnetic field and the 
two-photon state is still an eigenstate of the total linear 
momentum. However, in Sec. Ill we shall introduce an 
effective reduced density matrix poff both for the one- 
photon and two-photon states. 



In order to perform an actual calculation, we have to 
give a rule how to obtain the matrix elements P^s' (k) in 
Eq. 1)23(1 . From a physical point of view they can be ex- 
pressed as suitable linear combinations of the probability 
amplitude that a photon with definite momentum k and 
linear polarization Q!(k), is transmitted across a linear 
polarizer with a given orientation p. The next Section 
will be devoted to the search of an explicit expression for 

Pss'{^)- 



III. PHYSICAL POLARIZATION STATES 

In this section we study how a single-photon state 
changes when the photon crosses an arbitrarily oriented 
polarizer. Then we use this information to build an ef- 
fective reduced density matrix. 

A. Classical Polarization States 

A lossless linear polarizer j23| is a planar device which 
can be characterized by two orthogonal vectors: the axis 
n and the orientation p. The first vector is orthogonal 
to the plane of the polarizer, while the second one lies in 
that plane: p • fi = 0. From now on, with the sentence "a 
polarizer (p, n)" we shall indicate a linear polarizer with 
orientation p and axis n. 

For the moment we consider only classical fields, later 
we shall introduce the corres pon ding quantum operators. 
Following Mandel and Wolf |22|, we consider a polarizer 
(p, z) where p = x cos (3+y sin /3 and x, y , z form a carte- 
sian frame. Let and denote the incident and the 
transmitted electric field, respectively. We assume that 
E^ and E"^ are plane waves propagating in the direction 
z'. Then we can write explicitly: 



E^ = Ei,±' + El,y\ 

E^ = eI,s^ + eI,y', 



(40) 
(41) 



where x', y', z' are three orthogonal unit vector. If 6 and 
(j) are the spherical coordinate of z' with respect to x, y, z 
(where z is the axis of the polarizer), then 

x' = xcos^cos^ -|- y cos6'sin(/) — zsin^, (42) 
y' = — xsin0 -I- ycos(^, (43) 
z' = xsin cos (/) -f y sin0 sin -I- z cos 0, (44) 

and 9 < tt/2. The action of the polarizer can be found 
by requiring that the polarization e^^*^ of the transmitted 
field, lies in the plane defined by the polarizer orienta- 
tion p and the propagation vector z' of the impinging 
field [231 . In other words, e^^^ can be written as a linear 
combination of p and z': e'^^^ = cip -I- C2z'. The two real 
coefficients ci and C2 can be found by imposing normal- 
ization: je^P'l — 1, and orthogonality: e^P^ • z' = 0. The 
final result is 



DiP)e 



p - z'(z' • p). 



(45) 
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where D{P) = [1 — (z' • p)^]^/^. For later purposes, it is 
useful to write explicitly e'^'-' in terms of its components: 



p • X 

W) 



p ■ y 

Dm 



(46) 



Finally, we require that the transmitted field equals the 
projection along e'-^'^ of the impinging one: 



(E^-e(P))e(P). 



(47) 



Eq. (|47|l defines completely the action of the polarizer 
on the field. Now substituting Eq. (|1T|| and Eq. in 
the left and right side of Eq. H47|) respectively, and by 
using Eq. (|40|l . we obtain 



(48) 



where we have represented the impinging and transmit- 
ted fields as the row vectors and E^ respectively, and 
the transmission Jones 'IG] matrix T of the polarizer is 



T 



(x'-p)(y'-p) 

— m — 



(49) 



(x'-p)(y'-p) 

— m — 



/ cos^ 9 cos^ P cos 9 sin /3 cos /3 

1 — 6 cos^ j3 1— sin^ cos^ /3 

. cos 9 sin /3 cos /3 sin^ 

\ l-sin^ecos^/J l-sin^ S cos^ /3 



where (3 = (3 — (f>. It is easy to check that T is a real 
symmetric projection operator, and therefore TT = T. 
Moreover for 9 — Q — cj), Eq. H49|) reduces to the standard 
literature result for an "orthogonal" polarizer (p, z') (see, 
e.g., ref. [H, Eq. (6.4-20)). 



B. Quantum Polarization States 

Now we translate this in quantum language. Since we 
have introduced the generalized Jones matrix T for a 
linear polarizer, it is convenient to define the creation 
operators a\^^ of a photon with momentum k and linear 
polarization a = 1, 2 as 



(« = 1,2). 



(50) 



Since k/|k| = z' and the photon momentum is not af- 
fected by the polarizer, we define, for the sake of simplic- 
ity, olx' = Oki, o,y' = ak2- Therefore, from now on we shall 
always omit the momentum dependence in the expression 
of the one-photon states, only the bracket symbol |-) will 
remind that we are dealing with truly QED states. 

Again, we follow Mandel and Wolf 22] and introduce 
the operator vector field amplitude such that 



ay'Y . 



(51) 



The vector field amplitude sP" behind the polarizer, is 
determined by using the transformation law 

= flj/Ti/^/x' + a^>Ti,y,y', {i = x, y), (52) 

where summation over repeated indices is understood. In 
explicit form Eq. H52|l reads: 



where we have defined: 



bp — CLx' 



p • x' 



p ■ f 

D{P) 



(53) 



(54) 



Reminding that D'^{(3) — 1 — {z' ■ p)^, it follows imme- 
diately from Eq. (|54|) that the bp operators satisfy the 
canonical commutation rules 



[bp:bl]=l. 



(55) 



In order to clarify the relation between the quantum 
mechanical relation Eq. 1)54(1 and the classical Jones ma- 
trix Eq. (|49|) . let us briefly summarize our procedure 
in a more formal way. Let e(p) e be a given real 
unit vector depending from a set of real variables shortly 
indicated with p. As before, let and a^ denote the 
operator vector amplitudes of a plane-wave field before 
and after the polarizer, respectively. Then we define the 
annihilation operator bp in the same way as in Eq. 
by imposing 

a^ . e(P) = a^ • e^f) ee bp. (56) 
It is straightforward to show that: 

[bp^bl] = e(f)e(f)[al,(aj)t] 



(57) 



where summation over repeated greek indices is un- 
derstood. Moreover, as in Eq. (|50l) we have defined 
= aj • e|^"', {a = 1,2), and Eq. Q has been used. 
Eq. lf57|) suggests the introduction of a second rank sym- 
metric tensor defined as 



TP = p(p)p(p) 



(58) 



such that TrTP = 1. Now let us define e'^^'Hn terms of p 
as in Eq. 



(59) 



where D'^{p) = {Saf3 — kakp)paP(}. Then it is easy to see 
that the tensor coincides with the Jones matrix T in 
Eq.gni)- 

From Eq. H55|l follows that bp is a genuine bosonic 
operator, therefore we can use it, for different values of 
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/3, to calculate the states of the field behind the polarizer. 
In particular we define a one-photon state \ip{P)) as 



p • ic' 

W) 



\x'ii')) 



p ■ y 

Dm 



(60) 



where we have introduced the linear polarization basis 
states defined as 



l2/'(^')> ^ 4lO)- 

Moreover we define 

1^(0)) ^ St|0) ^ |x(z)), 
|^(V2)) ^ 5t|0> EE |y(z)>, 



(61) 



(62) 



where z refers to the axis of the polarizer and not to the 
photon momentum k = z'|k|. In fact, in a more complete 
way, one should write, e.g., |k, x(z)) for \x{z)). Clearly, 
Eq. (inH) reduces to Eq. (jHJ when = = 0. The two 
states Eq. (|62|l have unit length but are not necessarily 
orthogonal; in general we have 



cos^ 9 cos a cos /3 + sin a sin P 



[(1 



cos^ a sin^ ' 



'■9) (l-cos2/3sin2 6')]i/2' 
(63) 

where a = a — (j), (3 = P— cfi. From Eq. (|63|l follows 
that {tp{a)\tp{P)) = when /3 = ± arctan(cos 6*) and 
a = 2(f) + TT — p. In this case, when assuming arbitrary 
f}, the corresponding spatial orientations p(a) and p(/3) 
are not orthogonal: 



p(a)-p(/?) = -cos(2/3-20). 



(64) 



unless f3 ~ (j> + TT / 4. More generally, since (j) is arbitrary, 
we put = and we see that {ip{(3 + ^)\i/j{f3)) — when 



TT — arccos 



sin /3 cos /? sin 6* 



1 — cos^ /3sin^ 0(1 + cos^ 



, (65) 



which clearly reduces to 7r/2 when = 0. Therefore we 
conclude that is not possible to find a common orthogonal 
polarization basis for all values of 9. 

Finally we can answer the question posed in the end of 
the previous section. With the machinery we have built 
we can calculate, for example, the probability amplitude 
that an impinging photon in the state \x'(z')), is found 
behind the polarizer in the state |a;(z)); the result is: 



(x(z)lx'(z')) 



D{0) 



(66) 



More generally, we can calculate the non-unitary trans- 
formation matrix W^(k) as 



W(k) 



{xiz)\x'{i')) (x(z)ly'(z')) 
(y(z)|:.'(z')) (2/(z)|2/'(z')) 



COS COS fc 



Cl-cos2 0sin2 0)1/2 (i_cos2 0sin2 0)1/2 



sm cos £ 



cos <p 



(l-sin2 0sm2 0)1/2 (l-sin2 0sin2 0)1/2 

(67) 

The knowledge of M^(k) allows us to calculate all tran- 
sition amplitudes between both linear and circular po- 
larization states, the latter being obtained by forming 
suitable linear combinations of the elements Wij. Note 
that is a real matrix because we have chosen to rep- 
resent a linear polarizer; it is possible to show that in 
the case of an elliptic polarizer, complex phase factors 
appear p3 |. 



C. Effective Reduced Density Matrix 

In Sec. II we have shown that it is not possible in 
general to build a 2 x 2 reduced density matrix when 
many modes of the field are involved. However, looking 
at Eq. H29|) one is tempted to define an average 2x2 
density matrix p as 



N 



(68) 



This seems a reasonable definition since in a real experi- 
ment, the detectors automatically take averages over the 
photon momenta. Unfortunately this choice works only 
if the measured quantities do not depend on the momen- 
tum. In fact if, following the same line of thinking as 
above, we define an average 2x2 matrix projector P as 



1 ^ 



we can calculate its mean value as 



(P)=Tr(pP). 



(69) 



(70) 



This definition coincide with the original one given by 
Eq. (|28() only in the special case P{n) = P{m) \/n,m — 
1^ . . . , N. This result is hardly surprising, in fact by com- 
paring Eq. H19|l with. Eq. (|68|l one can easily recognize 
that p — pi- 

Let us try now a different approach by exploiting the 
analogies between classical and quantum optics. For a 
classical quasi-monochromatic light wave propagating in 
the direction k, a 2x 2 polarization density matrix p^s can 
be defined in terms of the measured Stokes parameters 
s„(z = 0,...,3) as mill 



1 

Peff- 2 



So + S3 Si 



IS2 



Si + IS2 So - S3 



(71) 
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Clearly the procedure one adopts to measure the Stokes 
parameters affects the actual value of Peg. Therefore, 
in classical optics, the polarization density matrix peff is 
understood as the measured density matrix, and different 
measurement schemes will lead to different matrices. Can 
we do the same in the quantum regime? As a matter of 
fact, when we have a well defined experimental procedure 
to measure the Stokes parameters of a beam of light, it is 
irrelevant whether the beam contains 10'^'^ or 1 photon. 
Therefore we regard the definition of Peft as given by Eq. 
(|71|l . as a postulate valid both in the classical and in the 
quantum regime. 

The quantum theory of light gives us the rules to calcu- 
late the Stokes parameters for both the one-photon [23 | 
and the two-photon states. We consider here only the 
one-photon case in some detail since the two-photon one 
is completely analogous. For a given momentum k and a 
polarizer axis z, let — {\k,x{z)), |k, ?/(z))} be the lin- 
ear polarization basis defined by two orthogonal polarizer 
orientation p(0) = x, p(7r/2) = y introduced in the pre- 
vious subsection. When k || z, in such a one-photon basis 
it is possible to represent the "Stokes operators" Si 01 
restricted to Tik, by the corresponding Pauli matrices: 
Si\Tik = (Ji, (i = 0, . . . ,3), where, e.g.. 



This step completes our calculation. The presence of the 
normalization constant Z^ff should not be surprising, it 
simply amounts to a renormalization of the Stokes pa- 
rameters Si with respect to sq: Si — > Si/sQ. 

Without repeating all the calculations, we shall give 
in the next Section directly the formula corresponding to 
Eq. (|75|l for the two-photon case. 



IV. TWO-PHOTON SCATTERING AND THE 
BELL-CHSH INEQUALITIES 

A. Two-Photon Scattering Matrix 

We consider now the following experimental configura- 
tion. A two-photon source emits a pair of polarization- 
entangled photons and sends them through two scat- 
tering systems 5*^ and Sb located along the photon 
paths. Two linear polarizers Pa and Pb are put in 
front of two multi-mode detectors Da and Db which can 
record both the two singles count rates and the coinci- 
dence count rate. 

The two-photon initial state emitted by the source is 
the Bell state 



= i(|y(z))(x(z)|-Kz))(y(z)|) 



(72) 





i 



When k z, the Pauli matrices transform as ct.; 
W^aiW. The mean values Si = {Si) can be calculated by 
using the total scattering density matrix pf = \ipf){ipf\ 
as 



Si = TT{pfSi), (i = 0, ...,3), 



(73) 



where the state \ipf) is given in Eq. 1)22(1 . Looking at Eq. 
(|72|l it is clear that all we have to calculate are the mean 
values of the four operators P^t (f, t — x,y) defined as 



P,r= J2 |k,fT(z))(k,T(z)|. 



(74) 



Note that the off-diagonal operators Pq-t t) do not 
correspond to physical observables, therefore are not Her- 
mitian. Finally, after comparing Eqs. H71I73|I and Eq. 
(|74|l . we can write 



(75) 



where Zcs = l/{{Pxx) + (Pyy)) is a normalization con- 
stant which ensures Trpcff — 1- Explicitly we have 



(77) 



where C± are complex coefficients such that |C+p + 
|C_p = 1, the subscripts A and B identify the two pho- 
tons, and the superscripts x and y denote the linear po- 
larization state. Moreover 

IV'?^) = |kF,a(kF))i„, ia = x,y; F^A,B). (78) 

The state is an eigenstate of the total linear momen- 
tum operator Kab = (8)1b + 1a'XiKb: 



J^AB\lpt) = (kA + kB)|V'i)- 



(79) 



So, at this stage, it is still possible to describe the state 
\ipi) in terms of a 4x 4 "polarization part" density matrix. 

The state of the pair after the scattering has oc- 
curred, can be written as 

1^/) = E E q)!"^' "(k)>A|q, /?(q))s, (80) 



where 



^a/3(k, q) = C+S^^Sa + C-S Sg^, (81) 



and 



5^ = (k,a(k)|k^,C(k^))i„, i^^x.y), 
Spri = (q,/3(q)|ki3,77(kB))i„, {'q^x.y), 



(82) 



Peff 



{P'xx) {Pyx) 



{Pxy) {Pyy) 



are the scattering matrix elements 27j. Moreover, K and 
(76) Q denote the sets of all scattered modes for photon A and 
B respectively, and a = a(k), (3 = /3(q). By inspecting 
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Eq. H80|l it is easy to see that the state |V'/) is not longer 
an eigenstate of the hnear momentum. Now, by repeating 
the same procedure we have executed for the one-photon 
scattering case, we introduce the two-photon generahzed 
reduced density matrix as 



P2 



f = ^ w;(k,q)|0(k,q))(0(k,q)|, (83) 



where we have defined 

|</.(k,q)) = J_^ ^vl;„^(k,q)|k,a(k))^|q,/3(q))B, 
VC(k,q) 

(84) 

where K^^, Qdb represent the sets of the scattered 
modes detected by detectors Da and Db respectively. 
Moreover, we have defined 



C(k,q) = ^|vl/„^(k,q) 



a,/3 



^2 = 1/ E C(k,q), 



(85) 



u;(k,q) = Z2C(k,q). 

As in the one-photon case, the operation of tracing with 
respect to the detected modes has reduced the pure state 
\'4'f) to the statistical mixture . Now it is clear that we 
can introduce a set of N^Nb^s pure state density matri- 
ces p(k, q) (4 X 4) whose elements are 



PQ/3,a'/3' (k, q) = 



^„ff(k,q)^;,^,(k,q) 
C(k,q) 



(86) 



where Na = dim/C^)^, Nb = dimQ/j^. Each of these 
matrices "lives" in a 4-dimensional Hilbert space Tikq: 



where (k e -fCi)^, q G Qdb) and {a, (3 = x, y). 

Not surprisingly, we have found a result analogous the 
one-photon scattering case, that is that a unique reduced 
4x4 density matrix does not exist. However, by using 
the methods developed in Sec. Ill we can introduce an 
effective 4x4 reduced density matrix p2off- In analogy 
with the one-photon case, P2cff can be written in terms 
of the measured two-photon Stokes parameters (i, j = 
0, . . . , 3) m as: 



Then, by following the same line of reasoning of the one- 
photon case, one can realize that it is possible to write 



(P2cff 1 



(89) 



where 

Pa'a = E |k, Q;'(zA))(k, a(z^)|, {a,a' = x,y), 



E 



|q,/3'(zB))(q,/3(zB)|, {P,P' = x,y), 



. (90) 

and 2'2eff such that Trp2cff = 1- With za and z^ we 
have denoted the axes of the two polarizers located in 
the paths of the photons A and B respectively. 

The calculation of p2off we have sketched above, closely 
resembles the previous one for the one-photon case. 
There is, however, an important conceptual difference 
between the two cases, as emphasized in Ref. |2^. In 
fact, in the two-photon case p2cff cannot be determined 
by local measurement only (each beam separately), but it 
is necessary to make coincidence measurement in order to 
account for the (possible) entanglement between the two 
photons. However, it is well known that entanglement 
properties of a bipartite system depend on the dimension- 
ality of the underlying Hilbert space [13 . IT3| , therefore 
the "measurement-induced" reduction from 4 x Na x Nb 
to 4 dimensions, may change the observed properties of 
the system. The problem of the determination of the ef- 
fective dimensionality of the scattered pair state, is at 
present under investigation in our group [Tsf. 



B. The Bell-CHSH Inequality 

We have just shown that, when in a two-photon scat- 
tering process we have a multi-mode detection scheme, 
the polarization state of the photon pair is reduced to 
a statistical mixture. We want to study the violation of 
the Bell inequality in the Clauser, Horne, Shimony and 
Holt (CHSH) form Ei], for that mixture. As usual the 
Bell operator Schsh is defined as Il7| 



Hkq = span{|k,a(k))^® |q,/3(q))B}, (87) 



CHSH 



a • cr ® (b + b') • (T + a' • cr (g) (b - b') • (T, (91) 



where a, a', b and b' are unit vectors in M"^. Moreover, cr 
is a vector built with the three standard Pauli matrices 
fi , (72, (73, and the scalar product a • cr stand for the 2x2 
matrix X]f=i o.i'^i- Then the CHSH inequality is 



|Tr(/56cHSH)| < 2. 



(92) 



In order to calculate explicitly Eq. H92|l it is necessary to 
know p which, in turn, depends on the specific scattering 
process considered. However, in our case, we want to 
show that the polarization-entanglement of a photon pair 
is degraded just because of the multi-mode detection, 
independently from the details of the process; therefore 
we shall consider a very general shape for p. 

Let ^(k, q) > denote the probability of a given phys- 
ical realization of the process. Then we can write 



E«^(k, q)|*kq)(«'kq| 
k,q 



(93) 
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where Ivfkq) represent an arbitrary polarization entan- 
gled state for which the photons A and B have momenta 
k and q respectively. This means that each time a pair is 
scattered, both photons will impinge on the correspond- 
ing polarizers with different angles determined by their 
momenta. So, for our purposes is enough to investigate 
the angular dependence of the entanglement of a single 
emitted photon pair, when at least one of the two pho- 
tons impinges with an arbitrary angle on the correspond- 
ing polarizer. In order to keep our treatment as general 




FIG. 1: Scheme of the two-photon analyzer described in the 
text. 5 is a source of Bell-Schmidt states, Pa and Pb are 
linear polarizers, and Da and Db are photo-detectors. The 
symbol & denotes the coincidence recorder. 

as possible, instead of considering a particular process- 
dependent scattered state, we focus our attention on the 
entanglement properties of the complete set provided by 
the Bell-Schmidt states TS 



Then we associate to each Bell-Schmidt state a well de- 
fined photon momentum pair (k^, k^) and we show that 
for each pure state the optimal choices of a, a',b 

and b' depend on (kyi, k^) and therefore it is impossible 
to find a choice which is simultaneously optimal for all 
the states in the ensemble given in Eq. (|93() . 

In order to demonstrate this, let us consider the detec- 
tion coincidence scheme shown in Fig. 1. An idealized 
source S emits photon pairs in the Bell-Schmidt states 
1$;). Two linear polarizers Pa and Pb are inserted in 
the paths of the two photons and two detectors Da and 
Db are put behind them. While Pa is put perpendicular 
to the momentum k^ of the photon A, the axis z of Pb 
is such that z • ks/lksj = cos9. 

Aravind has shown that the choices = 1, = 1, 
by = Q = by are optimal for all the therefore we 

make the same assumptions. The remaining components 
of the two vectors b and b' can be related to the physical 



orientations of the polarizer by writing, 

f = Tr(rcr) 

= {Tr(rai), 0, Tr(Ta3)}, 



(95) 



where f = b, b'. T is the polarizer Jones matrix as given 
in Eq. H49|l and Tr(rcr2) — because of the symmetry 
of T. Then we parameterize b and b' by introducing the 
two angles /3 and S respectively, in the Eq. (|95|l obtaining 



2 cos 9 sin /3 cos /3 cos^ 9 cos^ /3 — sin^ (3 



b' = 



1 — cos^ f3 sin^ 9 

2 cos 9 sin S cos S 



1 — cos^ (3 sin^ ( 
cos^ 9 cos^ S — sin^ 



1 - cos^ S sin^ 6* ' ' 1 - cos^ S sin^ 9 

' (96) 

where /? stands for /3 — and (5 for (5 — 4>. Now for each 
of the Bell-Schmidt states Eq. H94|l we choose the values 
for fi and b in order to maximize the violation of the 
Bell-CHSH inequality for 6* = 0, and calculate 



B,{9) = Tr(/3,6cHSH) 



(z = 1, 



-4), 



(97) 



where we have defined pi — |<I>i)(<I>i|. After a straight- 
forward calculations one finds that B\{9) — B2{9) and 
B3{9) = Bi{9), where 







\xa,xb) 


+ \yA,yB)), 






$2) 




\xa,xb) 


- lyA^ys)), 




B2{9) 








(94) 




$3) 




\xA,yB) 


+ \yA,XB)), 




B4{9) 


$4) 


= 7l( 


\xA,yB) 


- \yA,XB))- 







-4cos^ 



9\ 1 - (3 -I- 2 V2) cos 61 



2 J 1 (3 + 2V2) cos2 61 ' 



4cosM-^ l-(3-2V2)cos^ 



(98) 



2/ 1 (3 - 2V2) cos2 ( 



These functions are plotted in Fig. 2. It is clear that the 
optimal choice [Pq^^ , Sq^^) at = 0, is no longer valid 
when 9 increases and the degree of entanglement of skew 
photons appears to be reduced. However, one must real- 
ize that this loss of entanglement is an artefact due to our 
mismatched polarization detector. This means that it is 
still possible to find optimal values for {/3,5), but they 
will differ from the initial ones {9 = case). In order 
to show this explicitly, we have investigated the dynam- 
ics of points (/3°P*, i5°P') in the plane {(3,3), for varying 
9. The results are shown in Fig. 3 in the case of |<I'4): 
for the other states the results are qualitatively similar. 
When 9 increases passing from zero to 7r/2 , the points 
(/3°P*, (5°P*) move monotonically away from the central 
point (7r/2, 7r/2) along the line 5 = tt ~ (3 with different 
rates. Once S is fixed to the optimal value (5°p* = it — (3, 
one can follow the motion of as a function of 9. The 
dynamics is very simple and it is shown in Fig. 4. We 
list below the four functions P°^^{9), {i = 1,...,4) for 
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reference: 
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7T/8 7T/4 3/1/8 7T/2 
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FIG. 2: Average Bell-CHSH operator (Schsh) as a function 
of the polarizer "tilting" angle 8 = arccos(z • ks), for the 
Bell-Schmidt states |"I>4) and |'1?2)- The unit vectors b and b' 
have been chosen such that for both states (Bchsh) ~ 2\/2 
for 6 = 0. When 9 increases, the violation of the Bell-CHSH 
inequality decreases. The same curves have been obtained for 
|<l>i) (alike |<1>2)) and |$3), (alike 1$4)). 



to 




1 . 4 



0.0 



FIG. 3: The square boxes represent points in the plane (/3, 5) 
such that (Schsh) = ±2\/2 (+ lower series, — upper 
series) when the average is calculated with respect to |$4) 
for several values of 9. When 9 increases passing from zero 
(light grey) to about 7r/2 (dark grey), the "maximal violation 
points" move monotonically from light to dark along the line 

S = -K-P. 



f3°P\9) = -arctan[(l-\/2cos6')] , 
IT + arctan [(1 - V^cosO)] , 
arctan [(1 + \/2cos6')] , 
P2^\e) ^ n- arctan [(1 + \/2cos6')] . 



(99) 



Despite their simplicity, Fig. 4 and Eqs. tell us 

something important. We remind that the idealized ex- 
perimental scheme we have considered in Fig. 1 was in- 
troduced to study the behavior of an entangled photon- 
pair in a statistical mixture in which each photon-pair 
has a well defined momentum. While in the above anal- 
ysis we have considered as a free parameter represent- 
ing the polarizer axis, in a real scattering experiment 
is the angle at which one of the photons, belonging to 
the entangled pair, impinges on the polarization detec- 
tor. Then, each time a pair is scattered the two photons 
A and B will hit the detectors with arbitrary angles 9a 
and 9b respectively and the optimal polarizer orientation 
/J°P* will be different for each couple of angles. Therefore 
it is clear that we cannot simultaneously optimize /3 for 
all angles and the measured average degree of entangle- 
ment will be reduced independently from the scattering 
process considered. This completes our proof. Then we 
conclude that a conventional experimental setup for the 
measurement of the Bell-CHSH inequality, may fail to 
give the correct value for (yScHsn) when the measured 
state is a multi-mode scattered state. 



V. SUMMARY AND OUTLOOK 

The present paper aims to establish a theoretical back- 
ground for a future study of scattering processes by both 
chaotic optical devices 29] and random media. The main 
concern of this paper has been to demonstrate that in a 




(b) 



71/8 7t/4 3/T/8 71/: 



FIG. 4: (a-b) Dependence of the optimal value of (3 (such 
that (ZSchsh) = by the tihing angle 6. For 61 > 7r/4 

different values of 6 may require quite different values of /3. 
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scattering process the measured degree of polarization- 
entanglement of a photon pair is unavoidably decreased 
because of the multi-mode detection. At the core of this 
loss of entanglement process resides the correlation be- 
tween the momentum degrees of freedom and the polar- 
ization ones for the states of the electromagnetic field. 
In order to clarify the meaning of this correlation (or en- 
tanglement I ^6 have developed early on the paper a 
proper notation for the representation of the one-photon 
states of the electromagnetic field; this notation plays an 
important role throughout the paper. This introductory 
part also serves as a basis to show which dangers may be 
hidden behind the use of a misleading notation. In par- 
ticular we show that the use of a reduced density matrix 
obtained by blindly tracing out the momentum degrees 
of freedom, can lead to a wrong result when applied to 
the calculation of a polarization-dependent observable. 

The central part of this paper comprises two separate 
topics. The first topic consists in a careful analysis of 
the one-photon scattering processes. It is shown that a 
unique 2x2 reduced density matrix is an useless con- 
cept for the analysis of a multi-mode scattering process, 
and that more information than this is required. The 
second topic is how to build, within the QED context, 
the one-photon states selected by an arbitrarily oriented 
linear polarizer. The knowledge of these states allows 
us to introduce the concept of the effective reduced den- 
sity matrix which must be understood as the measured 
density matrix. 

The last part of this paper is devoted to a brief in- 



troduction to the subject of the two-photon scattering 
processes and to the investigation of the Bell-CHSH in- 
equality when, in a standard measurement setup, a po- 
larization analyzer is arbitrarily tilted. The violation 
of the Bell-CHSH inequality is explicitly calculated for 
the complete Bell-Schmidt set of polarization-entangled 
states. We show that, when in a two-photon scatter- 
ing experiment, the observer is ignorant about the mo- 
mentum distribution of the scattered photons, he cannot 
find an optimal orientation for the polarizers in order 
to maximize the measured violation of the Bell-CHSH 
inequality. However this does not mean that a scatter- 
ing process necessarily spoils the degree of entanglement 
of a given state, but instead just makes it not measur- 
able with a standard measurement setup. This naturally 
raises a question about the physical meaning of a com- 
putable degree of entanglement which does not coincide 
with the measurable one. This topic is currently under 
investigation in our group jl5| . 
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